Fate of 2D Kinetic Ising Ferromagnets and Critical Percolation Crossing Probabilities 
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We present evidence for a deep connection between the zero-temperature coarsening of the two- 
dimensional kinetic Ising model (KIM) and critical continuum percolation. In addition to reaching 
the ground state, the KIM can also fall into a variety of topologically distinct met ast able stripe 
states. The probability to reach a stripe state that winds a times horizontally and b times vertically 
on a square lattice with periodic boundary conditions equals the corresponding exactly- solved critical 
percolation crossing probability V a ,b for a spanning path with winding numbers a and b. 
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PACS numbers: 64.60.My, 05.40.-a, 05.50.+q, 75.40.Gb 

When an Ising ferromagnet with non-conserved spin 
flip dynamics is quenched from above the critical tem- 
perature to zero temperature, a beautiful coarsening do- 
main mosaic emerges [H-Q (Fig. 1). For finite systems, 
this coarsening of the kinetic Ising model (KIM) ends 
when the typical domain length reaches the linear di- 
mension of the system. What is the resulting final state? 
A naive expectation is that the ground state is ultimately 
reached because each microscopic spin update either de- 
creases or maintains the energy of the system. However, 
it is now appreciated that this lowest-energy state is not 
necessarily the final outcome. There exist a plethora of 
met ast able states, such as straight stripes in two dimen- 
sions [4j-[2( and more bizarre gyroid or "plumber's night- 
mare" states in three dimensions [7j, which are infinitely 
long lived at zero temperature. Once the system falls 
into such a state, the only escape route is via energy- 
raising spin flips. Since such events do not occur at zero 
temperature, there is no escape to the ground state. 

In the intermediate-time regime, where the typical do- 
main size substantially exceeds the lattice spacing but 
is much smaller than the system size, the domain mo- 
saic visually resembles the cluster geometry of contin- 
uum percolation [8]. This correspondence has sparked 
recent work on possible connections between these seem- 
ingly disparate models [1, In two dimensions, con- 
tinuum percolation is critical when the concentrations of 
both phases are equal This duality explains why 

the ground state corresponding to the majority phase is 
always reached (in the thermodynamic limit) whenever 
the initial magnetization is non-zero [4]. In this case, 
the majority phase percolates in all directions and in- 
evitably engulfs the entire system. The most interesting 
case of quenching from above the critical temperature 
corresponds to zero initial magnetization, so that the 
system in the intermediate-time regime is at the critical 
point of two-dimensional continuum percolation. 

The connection to critical percolation is extraordinar- 
ily fruitful because it allows us to understand why the 
system may fall into stripe states rather than ground 
states and it also predicts the probabilities of various 
outcomes [8[. For example, the probability to reach a 
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FIG. 1: Snapshots of coarsening in the nearest-neighbor ki- 
netic Ising model on a 1024 x 1024 square lattice with periodic 
boundary conditions at: (a,e) t = 200, (b,f) 1000, (c,g) 5000, 
and (d,h) 50000 after a quench from T = oo to T = 0. Top: 
evolution to (1,1) stripes (probability ~ 0.04); bottom: evo- 
lution to (2, 1) stripes (probability « 0.00015). 



state with vertical stripes (in the final state these stripes 
are straight and have widths > 2) equals the spanning 
probability Vo,i to have a path that spans the system in 
the vertical direction at the percolation threshold (and no 
spanning paths in other directions). The spanning prob- 
abilities Pq 7 i and Pi,o are exactly known |lll-ll3j. and 
this led to the prediction that the probability to reach a 
stripe state equals 0.3390 ... for the square with periodic 
boundary conditions, in agreement with numerical simu- 
lations [8] . (For free boundary conditions this probability 
^1-27 -Q.3558...). 
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Here we argue that the connection to percolation is 
much deeper and applies to a large family of positive- 
energy metastable states, of which straight stripes are 
merely the simplest members. We will apply the connec- 
tion to percolation to determine the probabilities to reach 
general stripe states that wind a times in one Cartesian 
direction and b times in the orthogonal direction for the 
two-dimensional KIM with periodic boundary conditions. 

To reveal this connection to percolation, it is essen- 
tial to extend the KIM to more distant interactions. On 
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the basis of universality [14j, the dynamical behavior 
should not depend on the interaction range in a funda- 
mental way. However, there are subtle but important 
interaction-range dependent effects that help expose the 
parallelism between coarsening and critical percolation. 
For convenience, we define range in terms of the Man- 
hattan metric, in which the distance between (0,0) and 
(x, y) is + By extending the KIM to weaker second- 
neighbor ferromagnetic interactions, one sees that the 
regular [1, 1]°° staircase shown in Fig.[2]become infinitely 
long lived. That is, there is an energy cost to flip any spin 
on either side of this staircase. The stability of this diago- 
nal staircase causes a stripe state that winds once around 
the periodic square (torus) in both the x- and ^/-directions 
to be infinitely long-lived at zero temperature. Similarly, 
extending the interaction range to third neighbors addi- 
tionally causes the [2, 1]°° and [1, 2] 00 stripe states to be- 
come infinitely long-lived (Fig. 1(h)). As the interaction 
range becomes infinite |15j , stripe states with arbitrary 
integer winding numbers (a, b) are infinitely long-lived in 
a square system. 
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FIG. 2: (a) [1, 1]°° staircase interface. With nearest-neighbor 
interactions (dashed square) interfacial spins can flip freely, 
but are stable with longer-range interactions, (b) [2, 1]°° stair- 
case. Interfacial spins can flip freely with (Manhattan) first- 
and second-neighbor interactions, but are stable with longer- 
range interactions. 



To make the quantitative correspondence between 
coarsening and percolation, we need to make use of sev- 
eral exact results for the probabilities for spanning clus- 
ters with specified top ologies on the torus at the per- 
colation threshold 13-13, 16-22|. As above, we label 
spanning clusters by their horizontal and vertical winding 
numbers, a and b respectively. Unique classes of span- 
ning clusters arise for each pair of values a, b ^ in which 
a and b are co-prime (i.e., a and b have no common di- 
visors). Stripes that are characterized by (a, b) and by 
(—a, —b) are equivalent and we therefore set a > 0. 

Let V a ,b(r) be the probability for a spanning cluster 
in continuum percolation with winding numbers (a, b) on 
a rectangle with periodic boundary conditions and with 
aspect ratio r = L y /L x . Here L x and L y are the linear 
dimensions of the system in the x- and ^-directions. This 
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We tacitly assume that r > 1; for r < 1, the spanning 
probabilities can be extracted from the obvious duality 
relation V a ,b(r) = Vb, a {\)- 

Let us consider the simplest crossing probabilities for 
a square L x L system: (i) Vo = Vq,i + Vi,o = 2Po,i? 
the probability for a vertical or horizontal stripe, (ii) 
Vi = Vi,i +V\-\ = 27^1,1, the probability for a stripe 
in the (1,1) or (1,-1) directions, and (hi) for n > 2, 
we define V n = 4P n ,i 5 the probability for a stripe in the 
4 distinct (±n, 1) and (±l,n) directions. The series in 
Eq. ([2]) converges rapidly in j and we also make use of 
the series representation of the Dedekind r\ function, 

W 2 )]- 2 = P-\1 + 2P 12 + 5P 24 + 10P 36 + ■■■), 
with p = e -7r / 6 , to give 

Jp 3 (i 



p-' + lp* 
2p 24 + Ap a 



lp 7 (l + 2p 12 ' "~ 2J ' «- 36 
fp- 



(3) 



32 „-lrr +3 ^ _|_ 2^12 i I i ,36 



5p 24 + 10p^ 



where the last line holds for all n > 2. These stripe 
probabilities are given to 4-digit accuracy in Table I, and 
our numerical data for the first three probabilities, which 
are accessible by simulations, are shown in Fig. [3] The 
agreement between theory and the simulation results is 
extremely good. 
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0.3388 


0.04196 


1.567 x 10" 4 


4.438 x 10" 9 


1.906 x 10" 15 



TABLE I: The probabilities V n for (n, 1) stripes on a square 
lattice for small n. 

A second natural set of interesting cases are diagonal 
stripes with tilt angle ±45° on an L x nL rectangle with 
periodic boundary conditions. Following the same cal- 
culational steps as those given previously for the square 
system, the series representation for the corresponding 
probability II n is given by 



n n 



8 ^7n 



[1 + 2p 



Yin 



2p 



24n 



(4) 



where again p = e ^' . From this expression, we nu- 
merically obtain the values shown in Table II (to 4-digit 
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accuracy). Our simulation data for U n for n = 2 and 
n = 3 are consistent with our predictions. For n > 4, 
II(n) is so small that is not practical to accurately mea- 
sure it by simulations. 
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FIG. 3: (a) Numerical values for the stripe probabilities Po, 
Pi , and P2 as a function of . The arrows indicate the exact 
values given in table I. 



Continuing this line of reasoning, the only stable stair- 
case in the (l,n) direction is n°°, where n = [l,n]. Sim- 
ilarly, (12) 00 is the stable staircase in the (2,3) direc- 
tion, (112)°° is the stable staircase in the (3,4) direc- 
tion, (122)°° is the stable staircase in the (3,5) direction, 
etc. The number of staircases going in the same direc- 
tion is infinite. For instance, the (1122) 00 staircase goes 
in the (2,3) direction, yet it is unstable. This instabil- 
ity indicates that there is another general rule to build 
allowed staircase interfaces [25j: only minimal represen- 
tations are stable. Analysis of stable staircases reveals an 
intriguing connection with the Farey sequences and the 
Stern-Brocot tree [24]. To illustrate it, we recall that for 
two neighbors in some Farey sequence, e.g. for \ and ^, 
their 'sum' is defined via the rule \ ®\ — = §, and 
this is taken as an indication that (23) 00 is the stable 
staircase in the (2,5) direction. 
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FIG. 4: Survival probability S(t) versus t on a semi- 
logarithmic scale for the KIM on a 64 x 64 torus with nearest- 



neighbor interactions (solid) and second-neighbor interactions 
TABLE II: II n for diagonal stripes on a L x nL rectangle. (dashed). 



An intriguing feature of arbitrary (a, b) stripe states 
is the intricate nature of the staircase interface between 
stripes when a and b are both large. The boundaries 
between the stripe states discussed thus far are either 
perfect straight lines (vertical and horizontal stripes) or 
a regular staircase that is inclined at 45° (see Fig. [2fa)). 
Stability with respect to single spin-flip dynamics im- 
poses severe restrictions on the form of these staircases. 
For example, a stripe with winding numbers (1,1) could 
hypothetically arise from a regular staircase that consists 
of alternating vertical and horizontal steps of length 2. 
However, such a staircase is unstable because the energy 
is decreased by flipping the corner spins. This length con- 
straint holds generally: adjacent vertical and horizontal 
segments in any stable staircase cannot both be longer 
than 2. Thus the only stable interface for (1,1) stripes 
is the regular staircase that we define as I 00 . This stair- 
case consists of the periodic sequence of building blocks 
1 = [1,1], in which [1, 1] denotes a unit-length horizontal 
segment followed by a unit-length vertical segment. 



The existence of an infinite variety of spanning paths in 
the KIM with infinite-range interactions also has intrigu- 
ing implications for the model with short-range interac- 
tions. Consider first the classic case of nearest-neighbor 
interactions. A useful diagnostic to detect metastable 
stripes with winding numbers a, b > 1 is to monitor 
the "survival probability" S(t), defined as the probabil- 
ity that there still exist flippable spins in the system at 
time t (the term flippable means that when such a spin 
is flipped, the energy of the system either decreases or 
remains constant). If there is a single coarsening time r 
that scales as L 2 , then one naturally expects that S(t) 
should asymptotically decay as e~ l l T . 

However, the actual behavior is markedly different 
(Fig. HJ , with the evolution of S(t) governed by two time 
scales [4]. The expected behavior where S(t) ~ e~ l l T 
holds until S(t) has decayed to approximately 0.05. At 
this point, the remaining configurations predominantly 
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have a (1,1) stripe topology (top line of Fig. 1). As indi- 
cated in Fig. [2fa) many of the spins along the interface 
that separates two diagonal stripes are in zero-energy en- 
vironments and can flip with no energy cost. The fluc- 
tuations of these freely-flippable spins lead to bulk diffu- 
sive motion for the interface. When two such diffusing 
interfaces meet, energy- lowering spin flips occur that ul- 
timately lead the system to the ground state. The decay 
of S(t) in this asymptotic regime is again exponential in 
time, but now with characteristic decay time that scales 
asL 3 0. 

For the KIM with (weaker) second-neighbor ferro- 
magnetic interactions, the vertical and horizontal stripe 
states, as well as the staircase [1,1]°°, are all stable 
at zero temperature in a square system. Thus at long 
times, any remaining metastable states are stripes with 
still higher winding numbers. This feature is reflected 
in the time dependence of S(t). The decay of S(t) in 
the second-neighbor KIM is qualitatively similar to that 
of the nearest-neighbor model, but the break in the de- 
cay now occurs when S(t) ~ 10 -4 (Fig. [4j). The long- 
lived states that remain beyond this break are predomi- 
nantly those with winding numbers (2,1) and (1,2) that 
ultimately relax to the ground state by the same inter- 
face diffusion mechanism as that for the (1,1) winding- 
number states with nearest-neighbor interactions. The 
probability for such stripe states to occur was previously 
found to be 1.567 x 10 -4 (Table I), consistent with the 
location in the break in the time dependence S(t). While 
these types of tilted stripe states are ephemeral (albeit 
with a lifetime that grows as L 3 ) when the interaction 
range is finite, they become permanent when the inter- 
action range becomes long-ranged. 

While we have focused exclusively on single spin-flip 
dynamics, the connection to critical percolation appears 
to be much more general. We have verified [25] that the 
same stripe probabilities arise in simulations of the time- 
dependent Ginzburg-Landau equation [l|-l3| as in the ki- 
netic Ising model. We therefore anticipate that the con- 
nection to percolation will be valid for any curvature- 
driven coarsening process with a non-conserved scalar 
order parameter. 

To summarize, we have presented evidence for a close 
connection between zero-temperature coarsening of the 
kinetic Ising model with arbitrary-range but decaying 
ferromagnetic interactions and critical percolation. The 
probabilities for the kinetic Ising system to evolve to 
a state that contains stripe paths with specified wind- 
ing numbers apparently coincides with the exactly-known 
spanning probabilities in two-dimensional critical perco- 
lation. This equivalence suggests that the domain ge- 
ometry of the kinetic Ising model coincides with that of 



continuum percolation at the critical point. 
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